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Abstract. The quark gluon plasma as produced in heavy ion collisions is exposed to early 
anisotropics in momentum space due to its rapid expansion. Such anisotropies can lead to 
non-abelian plasma instabilities, driven by unstable gluonic modes that can grow exponentially 
fast. These plasma instabilities can be simulated using a discretized version of gauge-covariant 
Boltzmann-Vlasov and Yang-Mills equations. In the stationary case, a turbulent cascade forms 
f> | in the strong-field regime, which is associated with an approximately linear growth of energy in 

collective fields. Early longitudinal expansion slows down the growth of unstable modes, and 
' the formation of soft gluonic fields depends crucially on the initial conditions assumed. 



1. Introduction 

Heavy ion colliders like RHIC or LHC create the quark-gluon plasma (QGP) in an anisotropic 
state. Due to its fast expansion, initially along the longitudinal but later also along the 
transverse directions, the plasma cools quickly and only exists for a duration of about a few 
tens of yoctoseconds (1 ys = 10~ 24 s). At early times and close to the center of the QGP, 
longitudinal expansion dominates over transverse expansion. This quickly leads to strong 
momentum anisotropies along the polar angle with respect to the collision axis pQ. Such early 
polar momentum space anisotropies could allow for a violation of the viscosity bound [2J, or lead 
to the emission of photon double pulses that are separated merely by yoctoseconds [3]. Most 
notably, polar momentum space anisotropies can induce Chromo-Weibel plasma instabilities, 
which are generalizations of Weibel or filamentation instabilities in ordinary electromagnetic 
plasmas It has been suggested early that these instabilities may play a fundamental role in 
the QGP [5j El [7]. In fact, already an infinitesimal amount of momentum space anisotropy causes 
the appearance of instabilities in collisionless plasmas [9]. In electromagnetic plasmas, the 
current filamentation instability develops magnetic islands on a fast electron time scale induced 
by the deflection of the original electron streams, which eventually leads to an isotropization 
of the electron distribution [10]. This fast isotropization on a time scale which is faster than 
ordinary perturbative scattering processes could explain the fast apparent isotropization and 
thermalization that is suggested by hydrodynamic modeling of the early QGP evolution [11[I12|. 

Non-abelian plasma instabilities can be studied using real-time lattice simulations |14] [15| [16] . 
The effective field theory for the collective phenomena at the soft scales is provided by gauge- 
covariant collisionless Boltzmann-Vlasov equations [T7]. The corresponding effective action is 
nonlocal and nonlinear [181 119] . but can be made local using auxiliary fields in the adjoint 




Figure 1. Comparison of average total 
field energy densities £ for SU(2) through 
SU(5) on linear scale in 3+1 dimensional 
simulations for a stationary system with 
anisotropy parameter £ = 10 [13J. 



Figure 2. The power spectrum for the 
electric distribution for various gauge 

groups SU(2) through SU(5) at late times 
80 < moot < 150. The distance between the 
lines is m OQ £\t ~ 11 [13J. 



representation W/s(x;v) which encode the fluctuations of the distribution function of colored 
hard particles |20j . They depend on a spatial unit vector which appears in the velocity 
v fi = p/y|p| of a hard particle with momentum p+ The Yang-Mills equations are given by 

D^A)F^=f, (1) 

where the current j u is calculated from the auxiliary fields 

The non-abelian Boltzmann-Vlasov equation for soft fields reduces to 

[vD(A)]Wf,(x;v)=F fh (A)v'Y, (3) 

with = — ig[A^, ■]. The scale of the hard particles drops out from these equations. An 
anisotropic distribution function /(p) is obtained by deforming an isotropic distribution /i so 
according to 

/(p) oc / iso (p 2 + frl) (4) 

with anisotropy parameter — 1 < £ < oo [8]. In order to solve these equations numerically, the 
3-dimensional configuration space is divided in a cubic lattice, on which a discretized version 
of the above non-abelian gauge-covariant Boltzmann-Vlasov equations is formulated |15j . The 
unit sphere of velocities for the auxiliary fields is described by a discretized set of unit vectors 
|21|, [T5] or by an expansion in terms of spherical harmonics |14j . 

2. Numerical results 

In the stationary case, the exponential growth of non-abelian plasma instabilities is limited in 
3+1 dimensions by non-abelian self-interactions. The exponential growth is limited by gluon 
self-interactions that are no longer negligible at a certain magnitude of soft fields. These self- 
interactions lead to a turbulence cascade which form a power-law distribution fk oc k~ u with 
a spectral index that turns out to be about u ~ 2. While simulations are usually based on 



the gauge group SU(2) [22^ [To] I23j . it has been confirmed that the same spectral index holds 
also in the QCD gauge group SU(3) as well as in higher gauge groups [13]. The systematics 
of the scaling with N c in the non-abelian regime is shown in Fig. [I] where the gauge groups 
SU(2) through SU(5) are compared. One observes that for different gauge groups the energy 
densities cease to grow at approximately the same value. In the following linear regime, larger 
gauge groups grow faster than smaller ones. Figure [2] shows the late-time behavior of spectra 
for various gauge groups. The spectra are multiplied by k 2 so that a scaling with u ~ 2 would 
correspond to a horizontal line. The slow growth at large momenta corresponds to the linear 
growth regime of Fig. [TJ 




Figure 3. Magnetic Weibel instabilities on 
an expanding background. The influence of 
different initial values is studied for modes 
with wave vector u = 30 in the direction of 
the anisotropy [23] . 



Figure 4. Chromoelectric, chromomagnetic, 
and total energy densities function of 

proper time. Proper time is normalized such 
that for Q s = 2 GeV each unit of Af 
corresponds to 1 fm/c [25] . 



A longitudinal expansion of the plasma at early times modifies the exponential growth, as 
there are two competing effects: On the one hand, plasma instabilities drive the exponential 
growth, on the other hand the longitudinal expansion suppresses the growth. The net effect 
is a reduction from a growth exponential in time to exponential in the square root of proper 
time. This has been numerically observed using the color glass condensate scheme [26} [27] and 
the discretized hard loop scheme [281 29] . In those simulations, an uncomfortable delay of the 
onset of plasma instabilities has been observed [29] . Collective fields decay and suppress Weibel 
instabilities at early times of the expanding plasma. It has been pointed out though that this 
suppression depends strongly on the initial conditions assumed [24] . In Fig. [3j the delay of the 
growth is displayed for various possibilities of initial conditions. Early work concentrated on a 
seed electric field with only IP (to) / [25] or a seed magnetic field with only A 1 {tq) ^ [2U] . 
Seed magnetic fields instead of seed electric fields increase the delay of plasma instabilities, as 
do mixed initial conditions for the fields. Surprisingly, if one considers initial fluctuations in the 
currents with only ^ 0, the delay is very strongly reduced [2"4] . 

This behavior has been confirmed in numerical simulations in 3+1 dimensions. Figure [3] 
shows various energy densities as a function of proper time for a longitudinally expanding 
system. Initially, the soft fields are depleted by the longitudinal expansion. Only after some 
time, the unstable modes overcome the depletion and all field components reach approximately 
the same magnitude. The growth rate is moderately reduced and transverse chromoelectric and 



chromomagnetic fields begin to dominate the energy density. Contrary to the fixed-anisotropy 
simulations, a saturation of the roughly exponential growth is not observed |25| . 

Since plasma instabilities may cause the isotropization in plasmas, the question remains how 
to measure the isotropization time. Direct photons would be a good indicator because they 
leave the QGP likely without further interaction. Under certain conditions, that is non-central 
collisions and photon production in a direction close to forward direction, non-trivial photon 
pulse shapes may be expected due to intermediate non-isotropic photon emission [3]. In extreme 
cases, these pulse envelopes may assume the shape of double pulses at the yoctosecond time scale. 
While it will not be possible to resolve such time structures directly [3U], it may be possible to 
observe the effect of such modifications to the pulse envelope through Hanbury Brown-Twiss 
|31] correlation measurements. A photon detector in the required forward direction may be 
installed during the ALICE detector upgrade by 2018 when the proposed Forward Calorimeter 
may be installed. With a few hundred photon pairs expected per year, such a measurement may 
be challenging, but not impossible. Thus, photons could provide valuable information about 
the earliest times of the plasma evolution, where gluon dynamics may be subject to plasma 
instabilities due to the rapid expansion along the beam axis. 

Concluding, one can state that for heavy-ion collisions at RHIC, non-abelian plasma 
instabilities probably may not have enough time to develop as they compete against the fast 
longitudinal expansion, but depending on the initial conditions, non-abelian plasma instabilities 
may play an important role at LHC energies [24, 25]. 
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